Invariant form of spin-transfer switching condition 
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We derive an invariant form of the current-induced switching condition in spin-transfer devices and 
show that for energy minima and maxima the "switching ability" of the current is determined by the 
spin torque divergence. In contrast, energy saddle points are normally stabilized by current-induced 
merging with other equilibria. Our approach provides new predictions for several experimental 
setups and shows the limitations of some frequently used approximations. 



High density electric currents induce magnetization 
motion and switching in nano-size metallic wires contain- 
ing alternating ferromagnetic and non-magnetic layers 
(Fig.[T]Ja)). This phenomenon is finding important appli- 
cations in computer memory and logic devices. Switch- 
ing is caused by the spin-transfer torque Tst^''^ which 
depends on the current, spin polarization, material pa- 
rameters and the geometry of the device. Once Tst is 
found, the magnetization dynamics can be obtained from 
the Landau-Lifshitz-Gilbert (LLG) equation. A sim- 
ple but often sufficiently accurate approximation is the 
macrospin model that assumes uniform magnetization of 
the layer M(r, = Mn{t), where M is the saturation 
magnetization value and n is a unit vector. In this case 
the LLG equation reads 
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where £(n) — {j/A'I)E{n), E{n) is the magnetic energy 
per unit volume, 7 is the gyromagnetic ratio, and a is the 
Gilbert damping constant. The spin-transfer torque is 
proportional to electric current, Tst ^ I- At / = vector 
n assumes an equilibrium position at a minimum 
of magnetic energy. A nonzero current has two effects: 
First, the spin torque gradually shifts the equilibrium 
away from its original position r\.eq{I = 0) ^ neg(/). 
Second, a stable equilibrium may abruptly turn unstable 
at a critical current /c, causing magnetic switching.^ 

Computation of Ic for a given equilibrium is a straight- 
forward though cumbersome mathematical procedure. It 
would be much simplified if one could find a single quan- 
tity that determines the "switching ability" of spin torque 
at a given equilibrium n^q . Often, it is implicitly assumed 
that the magnitude \Tst{'a.eq)\ itself is the relevant quan- 
tity, and, in particular, a sharp difference should exist 
between the destabilization of "coUinear" (rsi(neg) = 0) 
and "non-coUinear" {Tst{^eq) 7^ 0) equilibria. This view 
is certainly oversimplified, and it was argued^ that the 
critical current should also depend on the derivatives of 
Tst ■ Here we show that the "switching ability" can be in- 
deed introduced for extremum (minimum or maximum) 
energy points and in some cases for energy saddle points, 
and find explicit expressions for it. We apply our ap- 
proach to a number of experimental devices, obtain a new 
qualitative understanding of their dynamics, and clarify 



the limitations of some approximations. 

The LLG equation ([1]) can be equivalently written as 



(1 -I- a^) n F(n) = x(n) + an x T(n) , 
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with T = Tc + Tst, Tc = —[{de/dn) xn]. The equilibrium 
magnetization orientations lieq satisfy x(neg) = 0. Their 
stability can be investigated by linearizing the equation of 
motion. In spherical coordinates (0, 9) one decomposes 
F = F'^ecf, + F^eg in terms of the unit vectors e^, eg 
along the coordinate lines and obtains: 



59 



1 dF" _J_dF^\ /r, 
sine dd> sine dO \ "0 



dF" 



dF" 



69 



= D 



59 



(3) 



Stability of an equilibrium requires both eigenvalues of 
the "dynamic matrix" D to have negative real parts. This 
is equivalent to 



TrDiueq) < 0, detD(neg) > 0. 



(4) 



Noticing that matrix D is not covariant, we are led to 
introduce the related matrix of covariant derivatives 
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The trace of Dcov is an invariant quantity equal to 
Tri)...=divF. — (-(sin.FV^). 
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FIG. 1: (a) Typical spin-transfer device, (b) Angular depen- 
dence of the divergence and magnitude (inset) of the Slon- 
czewski's spin torque term^ for different polarizations P . 
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Crucially, D = Dcov at equilibrium points. Hence, an 
invariant condition divF < can be used instead of the 
first inequality in (j4|). The latter is relevant for treat- 
ing the equilibria corresponding to the energy extrema 
(minima and maxima). In this case det-D/=o > 0, with 
Tr_D/=o < at the minimum and TrD/=o > at the 
maximum points. Thus the condition of destabilization 
or stabilization is the change of sign of the trace, i.e., of 
divF. Using the relation between F and r and notation 
divF = V • F, where V operates in the n-space, we find: 

V • F = V • T - a[V X r] • n, (5) 

where [V x r] • n = —{dre/dcf) — d{s\i\9T^) / dd) / sin0. 
The general expression for the spin-torque created by a 
polarizer pointing along the unit vector s reads: 

Tst(ji,I) = uji g{n • s) [n X (s X n)] = ujjist{n) , (6) 

where uji = {j /M){hl /2eV), V is the magnetic layer vol- 
ume, e is the electron charge, and ^(n- s) is the efficiency 
factor;^ Using V • Tc = and [V x Tst] • n = we get 

V -F = V -Tst-aV XTa^ujiV -ist- aV^e . (7) 

Here the first term is proportional to the current and 
can lead to the sign change of the whole expression. We 
see that the switching ability is determined by the di- 
vergence V • fst that characterizes the angular depen- 
dence of the spin torque. Note that with Eq. ([7]) con- 
dition divF = can be viewed as a limiting case of 
the condition for the existence of a precession cycle (P), 
§{Tst ■ e±)dl — q;_^([V£ x n] • e±dl, where integrals 
are taken along the cycle^^ In terms of the bifurcation 
theory,^ local destabilization of the minimum points is 
the Hopf bifurcation which normally produces a stable 
precession cycle. 

Consider now the experimentally relevant case of small 
Gilbert damping, a ^ 1. Expression ([7]) shows that the 
critical current satisfies /c oc a and hence will also be 
small. Therefore at I — Ic the equilibrium point will be 
close to the zero current equilibrium, neq{Ic) — neq(O) -I- 
An with An oc /c oc a. Expanding ([7]) up to linear terms 
in OL we get an approximate stability condition 



with equality achieved at the critical current. Impor- 
tantly, all quantities in ([S]) are evaluated at the unper- 
turbed equilibrium point neg(O). In comparison, using 
conditions ([U one needs to perform an explicit calcula- 
tion of neg(/c) even in the case of the first order expansion 
in a (e.g. Ref. 00) ■ This welcome simplification stems 
from divF|/=o ~ ct holding for any n, while Tr_D|/=o ol 
holds only at neq(O). 

For Tst given by Eq. ([5]) one gets 



where is the angle between s and n. Representative 
graphs of V • fst(^?) are shown in Figdjb) for the Slon- 
czewski forn*^ of g{&). We observe that: (a) Divergence 
V • Tst can substantially differ from \Tst\^ i.e., the desta- 
bilization of noncoUinear equilibria may actually require 
larger current, (b) The switching ability V • fst vanishes 
at a critical angle (Fig[ljb)). Equation ((8]) predicts 
infinite critical current for the equilibrium points lying on 
the "critical circle" (CC) defined by ^((/), Q) = i?* (more 
precisely, at CC the approximation ([8]) breaks down and 
Ic is just large). The critical circle divides the unit sphere 
into two parts. Spin-transfer torque destabilizes the en- 
ergy extrema in one of them (which one - depends on 
the current direction), while in the other it makes them 
more stable, (c) The signs of /c are opposite for equilibria 
located on different sides of a CC. This circumstance is 
especially relevant when one considers different models of 
.g(z?). For example in the Slonczewski's case "i?* depends 
on the spin polarization P and varies from 'i?*(p=o) ^ ""/^ 
to 'i?*(p=i) = TT. In contrast, for a popular approxima- 
tion g = const, one has i?* = 7r/2 independently of P. 
The difference between the models becomes crucial for an 
equilibrium located between the respective CC's: a given 
current would have a stabilizing effect in one model, and 
destabilizing in another. 

As an illustration, consider a typical nanopillar'' with 
an [x, y) easy plane and an easy axis v J- z, so that 
e — ^ujp{nz)'^ — ^u!a{nv>)'^ — 7(H • n). Magnetic field is 
in-plane perpendicular, H ± i, H ± i^, and the polar- 
izer direction is s||i/ (Fig. [2]). In this setup the energy 
minima Mi and M2 are located at ±1/ at H = Q, move 
towards each other with increasing H , and finally merge 
with the saddle point L as the field reaches the easy axis 
anisotropy field Ha — The switching diagrams for 



(a) g^g(p) 2 (b) g = const ^ 
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FIG. 2: Critical circles and the "anomalous" stabilization re- 
gion. Upper panels: coUinear device with s||i/ and H _L z, 
H J- jy. The energy minimum points Mi and M2 move with 
increasing H as shown by the arrows. Critical circles shown 
for (a) generic (?(i?) and (b) g = const. Lower panels: switch- 
ing diagrams. In the regions Mi and M2 one equilibrium is 
stable, in B both are stable, and in P both are unstable. 
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a generic g{d) dependence (Fig. [2^) and the special case 
oi g = const (Fig. [^b) are qualitatively different^'^ with 
the former displaying the "anomalous" region (Fig. [2^). 
In our approach the "anomaly" is naturally explained by 
the fact that the minimum point Mi (H) crosses the crit- 
ical circle at H = iJ,. For g = const the minima never 
cross the critical circle, hence the anomalous region is 
absent. Note that the Tst{^) dependence produced by 
the g = const approximation is qualitatively similar to 
the actual one. Nevertheless, it does not lead to the cor- 
rect qualitative picture of switching when the equilibria 
of interest are close to the actual critical circle. 

A similar example is provided by the nanopillars with 
Hjli/ where magnetic field causes a crossing of the crit- 
ical circle by an energy maximum point. That crossing 
naturally explains the peculiar sign change of the corre- 
sponding critical current found in Ref . |3- 

Sensitivity to the g(9) angular dependence turns out 
to be of crucial importance for the interpretation of the 
"spin- flip transistor" (a nanopillar with s ± z, s ± u) 
precession experiment.^" Here the calculations with g — 
cons'^^'^^ forbid precession cycles (P) at zero magnetic 
field, but find them in external field H antiparallel to 
s. Based on this, Ref. [1^ interpreted the observation of 
precession at iJ = as an indication that an additional 
"field-like" term had to be introduced in Eq. ([5]) . Within 
the framework of our analysis, the absence of P states 
at zero field is due to the fact that zk. g — const the 
Af-points stay on the critical circle (Fig. [3^) and can- 
not be destabilized. The antiparallel field is required to 
shift the M-points away from CC. However, for general 
giiS), the A/-points are away from CC even at H ~ Q 
(Fig. [Sj^). They can be locally destabilized, producing P 
states by Hopf bifurcation without any field-like terms^iS 
The H = results of Refs. [H and[l3 are sensitive to the 
angular dependence 5(1?) in a manner that would be hard 
to foresee without the notion of a critical circle. 

To further demonstrate the power of the analysis based 
on Eqs. ^ and ([9]), consider experiments^ performed on 
a nanopillar device with an unusual "wavy" Tst (^3) depen- 
dence (Fig. 2^). In this case there are two critical circles, 
CCi and CC2, defined by the angles ^*i.2- At zero exter- 
nal field the energy minima Afi.2 fall into the regions of 



(a) g = const, H^O (b) g = g{^),H=0 




FIG. 3: Spin-flip transistor geometry, (a) For g — const local 
destabilization of the energy minima requires external field 
H to shift Afi,2 away from the critical circle CC. (b) For 
g = 3(1?) the minima are away from CC even at = 0. 
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FIG. 4: (a) "Wavy" r,i(i9) dependence (solid line, Ref. [TJ) 
and corresponding V ■ Tst (dashed line), (b) Critical circles 
CCi,2 and positions of minimum points Mi, 2 at an interme- 
diate value of field H. The energy maximum point is X. 



the same sign of V • Tst and can be destabilized simulta- 
neously, producing a precession cyclcj^ With increasing 
current, the cycle gradually approaches the energy max- 
imum point X. Eventually spin-transfer stabilizes that 
point-^ by closing the contour on it,''^^ The notion of crit- 
ical circles suggests an experiment capable of providing 
additional evidence for the "wavy" Tstii^) dependence. If 
an in-plane perpendicular H is applied (Fig.|4)D), the en- 
ergy minima Mi^2 are shifted towards the saddle point L. 
As "(5*1, 2 are not symmetric w.r.t. 7r/2, there will be an 
interval of fields where Mi had already crossed CCi and 
moved into the middle region, while Af2 remains in the 
left region. In this interval V • Tgt has opposite signs for 
A'/i,2 and normal switching between Mi and Af2 will be 
possible. Further increase of H will put both A/-points 
into the middle region, where they will be again destabi- 
lized by the same current direction However, in contrast 
with the H — case, now the same current direction 
will also destabilize X, so the P state evolution will be 
different ."^2, 

Let us now turn to stabilization of the saddle points. 
Here detZ)/=o < 0, so the process requires a change of 




FIG. 5: (a) Spin-flip transistor: s _L i, s ± 1^, H||s. The 
dashed lines show how the positions of the equilibria Mi and 
M2 change with increasing current and merge with the saddle 
L. (b) General in-plane directions of s and H. The saddle 
merges with one of the minima, while the other one asymptot- 
ically approaches s. (c-e) transformation of the fleld F during 
the merging of a saddle with two foci in case (a) . 
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(b) g = gW 




FIG. 6: "Magnetic fan" geometry, (a) At g = const the point 
M stays on CC until it collides with L, creating a large "OPP 
cycle" (dotted line), (b) At g = M is away from CC. 
Its local destabilization can create a small "IPP cycle" . 



sign of detD. An example is provided by a spin-flip tran- 
sistor where the spin torque attracts n to the saddle point 
and eventually stabilizes LJ^iii Notably, stabilization 
is always accompanied by a simultaneous discontinuous 
change in the nature of other equilibria. At H — the 
Af-points loose their stability just as L becomes stable. 
With the field parallel to s (note the difference with the 
antiparallel case discussed above) the current leads to a 
significant deviation of the Af -points from their initial 
positions (Fig. OJa)). At the critical current, Mi. 2 ap- 
proach L and merge with it, forming a stable center. We 
start by explaining why those simultaneous transforma- 
tions are not a coincidence. The saddle point is stabi- 
lized by becoming a stable center. As topological defects 
of the vector field F, saddles and centers differ in the 
winding number— which is a topological characteristic 
equal to n = — 1 for a saddle and n = 1 for a center 
or focus. Since the total winding number is conserved 
(the Poincare index theorem), a saddle point cannot be 
transformed into a center locally. The saddle-to-center 
transformation has to either proceed via merging with 
other defects (Fig. El^c-e)), or be accompanied by a si- 
multaneous change of nature of the far away equilibria. 

More insight comes from considering a generic case 
of H and s pointing in arbitrary in-plane directions 
(Fig. [3b)). Here L merges with one of the minima 
annihilating both equilibria, while the other minimum 
approaches s. Such merging is allowed by the winding 
number conservation and in fact the bifurcation theory^ 
shows that it is the most general case of the saddle-node 
bifurcation] i.e., the saddle point is normally not stabi- 
lized but rather destroyed in a collision with an energy 
extremum point. Stabilization happens only in special 



circumstances, such as s pointing exactly into L. In this 
case L remains an equilibrium for the arbitrarily large 
current and cannot disappear. That restriction produces 
the transcritical bifurcation^ where the energy extremum 
and the saddle exchange their nature in a collision. We 
find the critical current to be 
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7=0, Q=0 



5(0) 



(10) 



The above formula remains a good estimate in the case 
of a small misalignment between s and L. It also works 
for the spin flip transistor (Fig. [5^) even though here ad- 
ditional symmetries produce a more rare fork bifurcation. 

To sum up, the saddle point stabilization is associated 
with a merging or a close approach of equilibria, and thus 
can be detected without even calculating the dynamic 
matrix. A non-local bifurcation can only occur in devices 
of exceptionally high symmetry. 

Local destabilization and equilibrium merging are two 
alternative switching mechanisms which can compete 
with each other. Consider the s\\z "magnetic fan" 
experiments^ (Fig. [5]). Here the g = const approximation 
is special since the critical circle goes through the Ap- 
points. It predicts a merging of L and Mi, after which 
the system jumps into an "OPP cycle" (Fig. [B^). For 
angle-dependent g{i)) the CC is away from the Af-points, 
allowing for a competing scenario with local destabiliza- 
tion of Ml producing an "IPP cycle" (Fig. Eb)^ We 
estimate the critical currents as Lonipp) ~ aujp/ g' {-k /2) 
and oj I (OPP) ~ ^al gij^ 1"^)- Using Slonczewski's g and ex- 
perimental parametersi® we find loh^opp)/^i(ipp) ~ 0.2, 
which yields an OPP cycle scenario in accord with ex- 
periment. Furthermore, we predict that a device with a 
sufficiently large uja would manifest an IPP cycle. 

To conclude, we found that the ability of electric cur- 
rent to destabilize magnetization at its energy minimum 
(or to force it to be stable at its energy maximum) is de- 
termined by the spin torque divergence. We also showed 
that the sphere swept by the magnetization vector is di- 
vided into stabilization and destabilization regions by 
critical circles. Finally, we demonstrated that saddle 
points are stabilized via a topologically distinct route of 
merging with other equilibria, and discussed the compe- 
tition between such a merging and local destabilization. 
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